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I. INTRODUCTION 

Current semiconductor technology can produce a class of an electric modulation in a two-dimensional system. In 
a two-dimensional .electron system under a magnetic field such a modulation causes opite interesting structure in 
various observablesEl. It is also used to examine various aspects of quantum Hall systemcl. The system subjected to^f 
periodic potential which is described by finite number of harmonics has been extensively studied by many authorsutJ. 
The system subjected to a cosine potential is equivalent to the nearest-neighbor (NN) tight-binding model with a 
magnetic flux. Hofstadter computed its spectrum of butterfly-shape and discovered a multi-fractal structure!! Other 
periodic potential may also have a self-similar structure in a quantum region. Experimental observation of such a 
structure is a challenging theme. The Hall conductance is a significant observable to extract physical properties of 
such a system. Thouless et al. showed that in the system subjected, to a periodic potential the Hall conductance is 
written by the Chern number of the Bloch function in a gap regionen3. In the cosine potential case the Chern number 
is given by a solution of the Diophantine equation. The spectra and the Hall conductances in the next NN modelcl and 
in periodic potentials which are described by finite number of harmonicsd are also obtained. However, the problem 
of a periodic short-range potential remains to be solved. Recently, a remarkable progress wag -made. That is, its 
spectrum in the lowest Landau level was computed with sufficient accuracy to see self-similarityfl'E3. However, effects 
of higher Landau levels and Landau level mixing is not known well. In addition the Hall conductance of subbands is 
unknown. A suitable base function is necessary for this calculation. Many authors have used the Landau function 
or the eigenfunction of the angular momentum. The former one extends in one direction of spatial coordinates and 
localized in the other direction. This is suitable for the potential of one-dimensional translational invariance. The 
latter one is localized in the radial direction and has a rotational invariance. This is suitable for a rotational invariant 
potential. I— . .— . 

The magnetic von Neumann latticeEj'Eij is another basis in a magnetic field and can be constructed independently 
of the gauge. The basis consists of direct products of coherent states in the guiding center coordinate space and 
harmonic oscillator eigenstates in the relative coordinate space. The von Neumann lattice basis is a set of localized 
states. It turns into a set of extended states by the Fourier transformation. This basis has the following desirable 
properties: (i) A two-dimensional lattice translational invariance exists, (ii) Lattice structure of the von Neumann 
lattice can be changed according to a problem. (iii)A modular invariance of the von Neumann lattice exists. Owing 
to the first property, the Hilbert space of one-particle states is specified by the Landau level index and the lattice 
momentum. The second property is desirable in solving various periodic potential problems. The modular invariance 
is a key to develop a field theoretical formalism based on the von Neumann lattice. In our previous paper we showed 
advantages o£,the von Neumann lattice by studying the spectrum of a periodic short-range potential in the lowest 
Landau levellS. A field theoretical formalism is necessary when we study various interactions in a systematic way. In 
fact the field theoretical formalism based on-Uae.von Neumann lattice allow us to prove the exactness of the quantized 
Hall conductance by a topological formulatllta. The topological formula is written by the winding number of the 
full propagator in the momentum space. The weak localization, QED effects, and electron-phonon interaction do not 
alter the quantized Hall conductance. However, a periodic potential is not included in the proof. 

In the present paper we study dynamical properties of the system subjected to a periodic potential using the von 
Neumann lattice. We first construct explicitly three sets of one-particle states based on the von Neumann lattice. 
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That is, the set of coherent states, the set of momentum states and the set of Wannier states. The coherent state 
is localized on a lattice site and is not orthogonal. The momentum state is extended and orthogonal. Its analytical 
expression in the spatial coordinate space was not known. We give it-4n terms of the theta function for the first 
time. The-. Wannier state is localized on a lattice site and is orthogonalEJ. Recently, this state was discussed by two 
groupsll3il3. However, their results are restricted to a rectangular von Neumann lattice. We study the analytical 
form of the Wannier state in a general von Neumann lattice. Using the set of momentum states we develop a field 
theoretical formalism in a complete and strict manner. We pay attention to a peculiar property of the system in a 
magnetic field. In this system, we can not diagonalize both the current and the energy with respect to the Landau 
level index simultaneously. This imply that the usual relation Vi = J^- does not hold in this system. However, the 
gauge invariance allow us to prove a similar relation as the Ward-Takahashi identity. This identity is used to derive the 
topological formula of the Hall conductance. In a periodic potential system the spectrum consists of Landau subbands 
and the range of the lattice momentum becomes narrow according to the period of the potential. The propagator 
has additional index which reflects the subband structure. We prove that in this system the Hall conductance is 
written by the winding number of the propagator. Furthermore, this winding number can be rewritten by the Chern 
number of the eigenfunction. This Chern number is the same as the one obtained by Thouless et al.u. As an example 
we calculate the Chern number of the system subjected to a periodic short-range potential. The result gives some 
insights to the Wigner crystal phase of the quantum Hall system. 

The content of this paper is as follows. We construct three sets of one-particle states based on the von Neumann 
lattice in Sec. II. A field-theoretical formalism based on the set of momentum states is developed in Sec. III. The 
Ward-Takahashi identity and the topological formula of the Hall conductance written by the winding number of 
full propagator are shown. A periodic potential and defect are studied using the formalism in Sec. IV. The Ward- 
Takahashi identity and the topological formula are generalized to a system subjected to a periodic potentials in Sec V. 
The relation between the winding number and the Chern number are clarified. The Hall conductance of the system 
subjected to a periodic short-range potential are calculated with this formalism. Summary and discussion are given 
in Sec. VI. 



II. THE VON NEUMANN LATTICE BASIS 



In a two-dimensional system under a uniform perpendicular magnetic fielcLS, we introduce two sets of coordinates, 
i.e. the guiding center coordinates (X, Y) and the relative coordinates (£,77)13: 

t=(eB)- 1 (p y + eA v ), r? = -(eB)- l (p x + eA x ), (2.1) 
X = x-£, Y = y~ v , (2.2) 

where B — d x A y — d y A x and eB > 0. Each set of coordinates satisfies the canonical commutation relation 

and two sets of coordinates are commutative. Using these operators, a one-body Hamiltonian for a free charged 
particle is written in the form 

H a = \muj 2 c (e+V 2 ), (2.4) 

where uj c = eB/m. Hq is equivalent to the Hamiltonian of a harmonic oscillator and the eigenvalue is solved as 
follows: 

H \fl)=E l \f l ), ^=^ c (/+i), Z = 0,l,2,.... (2.5) 

This energy level is called the Landau level. Since Hq is independent of X and Y, the phase space (X, Y) corresponds 
to the degeneracy of the Landau level. It is convenient to use the coherent state defined by 

(X + iY)\a 

z mn = (muj x + nu y )a, 

where m,n are integers and uj x ,uj v are complex numbers which satisfy 
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Im[w> y ] = 1. 



(2.7) 



z mn is a point on the lattice site in the com plex plane; an area of the unit cell is a 2 . We call thisdattice the magnetic 
von Neumann lattice. Under the condition (2.7), the completeness of the set {|a mn )} is ensuredlia. The coherent state 
\&mn) is constructed as follows: 



|cw) = e -(" l +"+ m ")+^ t£ ^-^%)|aoo), 
A=^(X + i Y), [A,Ai] = l. 

The coherent states are not orthogonal, that is, 



(2.8) 



y-m-\-m' ,n-\-n' \^m' ,n' 



^7r(m+n+mn)- 



(2.9) 



Thus, a translational invariance exists. The Hilbert space of one-particle states is spanned by the state \fi (g) a mn }- 
The expectation value of the position of this state is a coordinate of the site of the magnetic von Neumann lattice 
and its mean square deviation is (I + l)a 2 /2ir. Hence, the state \ fi (£> a mn ) is a localized state. The wave function of 
the state \fi (g) a mn ) in the spatial coordinate space is given by 



(x|/i <g> Or, 



1 / 7P , Z - 



Ye 



r \z-Zmn \ 2 — iiT7n — +i7rn — +iA(x) 



(2.10) 



Here, £i = XjW^ 1 with the matrix W defined by 



W 



Re[w x ] Im^] 
Refo-iy] Imfwy] 



detW = 1. 



(2.11) 



In terms of x, y, z = x + iy — xLu x + yLu y . Thus, a site of the von Neumann lattice corresponds to integer x/a, y/a value. 
The function A(x) represents a gauge degree of freedom of the eigenfunction. For example A(x) = corresponds to 
the symmetric gauge. 



(a) 1=0 




(b) 1=1 




FIG. 1. The probability density a 2 \(x\fi ® aoo}| 2 with I — 0, 1 for the square von Neumann lattice. 

In Fig.l we show the probability density a 2 \(x\fi®aoo}\ 2 with I = 0, 1 for the square von Neumann lattice. Localization 
around x = is clearly seen in this figure. Since the probability density is a function of x — ma and y — na, the 
correspon dent figure of \ fi a mn ) is obtained by sliding the coordinate as (x, y) — > (x — am, y — an). 

As Eq. (2.9) has a translational invariance, an orthogonal basis can be obtained in the momentum space. Fourier 
transformed states denoted by 



Or 



(2.12) 



are orthogonal, that is, 



(ap|<v) = a(p)J2(^) 2 S(p - P' - 2ttN). 



(2.13) 



N 
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Here, N = (N x ,N y ) is a vector with integer values and p = (p x ,p y ) is a m omentum in the Brillouin zone (BZ), 
that is, \p x \, \Py\ < 7T- The function a(p) is the Fourier transform of Eq. (2.9) and calculated by using the Poisson 
resummation formula as follows: 



«(P) = Pfa)*P(p), 

P(p) = (2Imr)3 e l 3T p ^i( 



Px + rp y 
2ir 



r), 



(2.14) 
(2.15) 



where $i(z|r) is a theta function and the moduli of the von Neumann lattice is defined by r = ~lj x /lo v . The magnetic 
von Neumann lattice is parameterized by r. To indicate the dependence on r, we sometimes use a notation such as 
/3(p|r). For t = i, the von Neumann lattice becomes a square lattice. For r = e' i27r / 3 , it becomes a triangular lattice. 
Some properties of the above functions are presented in Appendix A. Whereas a(p) satisfies the periodic boundary 
condition, (3(p) obeys a nontrivial boundary condition 



(3(p + 2ttN) 



_ p i<p(p,N) 



P(p), 



(2.16) 



where <f>(j>, N) 
follows: 



tt(N x + N y ) — N y p x . We can define the orthogonal state which is normalized with (5-function as 



\Pp 



/3(p)e*x(p) ' 



(2.17) 



The function x(p) reflects an ambiguity with respect to the normalization. This ambiguity implies another gauge 
symmetry of this system. Later, we identify this symmetry as the gauge symmetry of constant magnetic field in BZ. 
It should be noted that the state |an) is a null state, that is, J2 m \a mn ) = dl3, because (3(0) = 0. Since \ao)/(3(0) 



is indeterminate in Eq. (2.17), \(3 ) is defined by limp^g \Pp)- The Bilbert space of one-particle states is also spanned 
by the state \fi ® (3 P ). We call the state \fi <S> (3 P ) the momentum state of the von Neumann lattice. The wave function 
of the state |// ® (3 P ) in the spatial coordinate space is given by 



(x|/, ® (3 P ) = J ^(f ) l (-2d z 



z) l \(3*(px 



2TT^,p y + 27T-)e lp " 
a a 



(2.18) 



The probability density |(x|/j (g) /3 P )| 2 is invariant under the translation (x,y) — > (x + aN x .y + aN y ) with integer 
N x ,N y . Thus, the momentum state is an extended state. 



(a) 1=0 




(b) 1=1 






FIG. 2. The probability density a 2 |(a:|/ ; ® (3 )\ 2 with 1 = 0,1 for the square von Neumann lattice. 

In Fig. 2, we show the probability density for the lowest two Landau levels. In Fig. 2 the von Neumann lattice is the 
square lattice and we choose p = 0. The probability density with a non-zero momentum p is obtained by sliding the 
coordinate as (x, y) — > (x + ap y /2n, y — ap x /2i:). 

Since \(3 p ) is orthogonal and normalized with 5-function, its inverse Fourier transformation leads to an orthonormal 
state. That is, the state defined by 



\Pn 



BZ 



d 2 p 
(27T) 2 



(2.19) 
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satisfies 



'I I 



— Swim/ 



mm u nn' • 



(2.20) 



The state \/3 mn ) is completely different from the coherent state \a mn ) due to the normalization factor /3(p)e lx ( p ' . The 
relation between the coherent state and the orthonormal state is given by 



\ft„ 



m — to ,n — n )\Q. r 



(2.21) 



where G(m,n) is the inverse Fourier transformation of 1/ f3(p)e lx< ^ p \ It has a long tail proportional to (to 2 +n 2 )~ 1 / 2 . 
The Hilbert space of one-particle states is also spanned by the state \fi (g) (3 mn ) which is an orthonormal localized 
state and has the center around z = z mn in the spatial coordinate space. However, it also has a long tail away from 
the center. We call the state \fi £g> /3 mrt ) the Wannier state of the von Neumann lattice. The wave function of the 
state \fi £g> fimn) in the spatial coordinate space depends strongly on the phase x(p)- A simple form is obtained if we 
choose x(p) — PxPy/Zft + (Py + 7r )/2- In this case the wave function becomes 



(x|/j ® P„ 



e iA(x) rr a n 
V7T (9- - 29 -* + ^ 

a V 11 Zir a z 

.(2I mT )3, /J_ - l 7rr*(f-m) 2 -27rim|+i7r(m+n+i)+ 4 7rf| C ^(fe-j 
' V ' V ir* J_„ 2lT 



+n+i+T*(f-m)) 2 



(2.22) 



The similar result was obtained by Zak for a rectangular von Neumann lattice in the lowest Landau leve!0. In Fig. 3 
we show the probability density a 2 (x|/; ® /3 00 ) with I = 0, 1 for the square von Neumann lattice. 



(a) 1=0 




(b) 1=1 




FIG. 3. The probability density a 2 \(x\fi <g) poo}\ 2 with Z = 0, 1 for the square von Neumann lattice. 



In Fig. 3 we show the probability density a 2 (x.\fi ® /3qo) with i = 0, 1 for the square von Neumann lattice. Asymmetry 
between x and y-direction is owing to the phase choice. The center of the wave function (x|/i (8) /?oo) is at [x, y) = 
(0, l/2)a in contrast to a naive expectation that it is at (x,y) = (0,0). This is not restricted to the case /?oq. Since 
the probability density |(x|/; ® (3 mn }\ 2 is a function of x — am, y — an, the slide of the cepter occurs for all (3 mn . We 
should note that the Wannier state found by Rashba et al. is slightly different from ouraia. They used |/3(p)| as the 
normalization factor. As a consequence, the probability density of the wave function (x|/o (8 /?oo) for the square von 
Neumann lattice has a center at z = and is rather symmetric_than ours. Their wave function behaves as 1/r 2 for 
large r which is the critical behavior in the Thouless's criteriontj. Our wave function behaves as 1/x for large x and 
is exponentially dumped for large y. This behavior is also the Thouless's critical behavior. 



III. FIELD THEORETICAL FORMALISM AND TOPOLOGICAL FORMULA OF HALL 

CONDUCTANCE 

In the preceding section we obtain three sets of one-particle states based on the von Neumann lattice, that is, 
the coherent state \fi ® a mn ), the momentum state \fi <St-ffp) and the Wannier state |// (g> (3 m n)- A field theoretical 
formalism based on the coherent state was developed inHl. A formalism based on the momentum state or on the 
Wannier state was partially used inc2l. Here we develop the field theoretical formalism based on the momentum state 
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in a complete and strict manner for later convenience. From now on, we denote |/z<8>/3 p ) as \l,p) and choose x(p) = 0- 
The reason to prefer oufnchoice (x(p) = 0) is simpleness of the form for the density operator p(k) discussed later. If 
we use the definition ok-3, p(k) becomes a complicated form. 
We expand the electron field operator in the form 

V>(x)= / — ^^^(pXx^p). (3.1) 

bi(p) satisfies the anti-commuting relation 

{bi(p),b\,(p')} = Si, v ^(2^) 2 J(p - p' - 2^)e^'' N \ (3.2) 

JV 

and the same boundary condition as /3(p). b\ and bi are creation and annihilation operators which operate on the 
many-body states. The free Hamiltonian is given by 

Ho = J dW( x) £ v,(x) = E / Bz (02^(P)MP)- (3-3) 
The density and current operators in the momentum space are 

p(k)=/ d 2 xe- jk - x V^(x)V(x), 

j(k) = J d 2 xe^ k - x (-^)V t (x)(^y^ + ieA(x))V»(x). (3.4) 
Using a modular transformation of /3(p), which is given in Appendix A, = (ft,j) becomes 

i M (k) = / ^E 6 Kp) 6 HP + «k)(/4K,e-^}|/ i ,)e^^( 2 ^)«. (3.5) 

Here, — (l,—u> c r),L) c £), and hi = Wijkj. The explicit form of (/; \e~ lk '^\fi>) is given in Appendix A. The phase 
factor in Eq. ( [3.5] ) can be written as exp(i fp +ak Ai(p)dpi), where the path is the straight line from p to p + ak and 
Ai(p) — (p y /2tt, 0). If we choose a non-zero x(p)> then Aj(p) beco mes A j — di\, i-e. a gauge transformation of Ai. 
Therefore, Ai(p) is a gauge field on BZ and the phase factor in Eq. (3J5) represents a unit flux on BZ. We note that 
the simple form of j'^(k) is due to our convention x(p) = 0- F° r later convenience, we derive the density operator in 
the spatial coordinate space. Let us define the polynomial /;;/ (fc) as 

(fl\e~ iM \fv) =fw(k)e- s & (3.6) 
The explicit form of fw (k) is given in Appendix A. The density operator can be written as 



a ijfJBZ (2tiT {2ttY idx 



x(3( Px - 2v V -,p y + 2n-)e-^ p y[3*(q x - 2^,q y + 27r-)e 1 ^. (3.7) 
a a a a 

Since p{x) = tp' (x)ip(x), we obtain a useful identity which relates the normalization factor /3(p) with the coordinate 
representation of the state |Z, p) as 

a%p\x)(x\l',<$ = fw(4r)P(P* - 2 *~,Py + ^-)e-^P*(q x - 2^,q y + 2^-)e 4 K (3.8) 
iox a a a a 



As a consequence of the above identity, we obtain Eq. (2. IS) for example. 

The free Hamiltonian 7i is diagonal in the above basis. However, the density operator is not diagonal with respect 
to the Landau level index. This basis, which we call the energy basis, is convenient to describe the energy spectrum 
of the system. In another basis, Ho is not diagonal and the density operator is diagonal. This basis, which we call the 
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current basis, is convenient to describe the Ward-Takahashi identity and the topological formula of Hall conductance. 
There is no basis in which both the Hamiltonian and the density are diagonal. This is one of peculiar features in a 
magnetic field. 

The current basis is constructed as follows. Using a unitary operator, we can diagonalize the density operator with 
respect to the Landau level index. We define the unitary operator 

Ul(p) = (fi\^ a -^™«\f v ). (3.9) 

By introducing a unitary transformed operator fe/(p) = Y],, Uw(j))bi(p), the density operator is written in the diagonal 
form and the current operator becomes a simple form: 

p(k) = / -|^^(p)Mp + ak). 
Jbz \ Z7r ) ; 

j(k)= / -^^^(pHv+^^p+^K^^P + ak). (3.10) 

bi and b\ satisfy the anti-commutation relation and boundary condition 

{Mp)^(p')} = E( 2 ^) 2(5 (P - P' " 27rN)A„,(N), 

N 

5j(p + 27rN) = 5^A W /(N)6i,(p), (3.11) 

v 

where A is defined by A(N) = (—l) N:c+N yU(2Tv'N) and is independent of p. However, it shuffles the Landau-level 
index. 

Here we review briefly the Ward-Takahashi identity and the topological formula of Hall conductance using the 
current basis. We consider various interactions perturbatively. It is convenient to define the one-particle irreducible 
vertex part r M as 

<j»{q)ll{p)lW)) = (2n) 3 5(p + Q-p>)SuApK ll2 (p,P + Q)S l2l ,(p + Q), (3.12) 
where Q M = (qo,aq x ,aq y ) = t^ v q v is a linear combination of q M and S is the full propagator defined by 

(bi(p)b\,( P ')) = (2Tr) 3 5(p-p')S w (p). (3.13) 

Between S and L M the Ward-Takahashi identity is satisfied. The identity has crucial roles in the following derivation 
of the topological formula of Hall conductance. The Ward-Takahashi identity in this case becomes 



dp" 



T li {p,p) = t\^^. (3.14) 



In a theory without a magnetic field, the Ward-Takahashi identity gives a relation that the state of the dispersion e(p) 
moves with the velocity S ^ > ) ■ However in a magnetic field, we can not diagonalize both the current and the energy 
simultaneously. Therefore, the Ward-Takahashi identity Eq. ( |3.14 ) does not imply the relation. 

The Hall conductance is the slope of the current correlation function Ti l - iv (q) at the origin and is written as 

<Txy = ^ p d p n^(q)\ q=0 . (3.15) 

If the derivative d p acts on the vertex with the external line attached, its contribution becomes zero owing to the 
epsilon tensor. Therefore, the case that the derivative |acts, on the bare propagator is survived. At this point it is 
proved that only diagrams of Fig. 4 do contribute to a- xy E3'EJ. 
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FIG. 4. Feynman diagrams which contribute to a xy . 

In Fig. 4 the dark disk of the vertex and the dark disk of the propagator are the one-particle irreducible vertex part 
and the full propagator. These satisfy the Ward-Takahashi identity. Thus, <j xy is written as a topologically invariant 
expression of the full propagator: 

cr xy = ^ d 3 pe^ p tr ^S- 1 (p)S(p)dJ- 1 (p)S(p)d p S- 1 (p)S(p)) (3.16) 



BZxS 1 



Here, the trace is taken over the Landau level index and the po integral is a contour integral on a closed path which 
is drawn in Fig. 5. 

t vPo 



FIG. 5. The contour line of po- Bold lines in the real po axe represent Landau bands. The Fermi energy Ef is located in the 
gap region. 

Thus, we denote S 1 as the integration range. _ 
If the Fermi energy lies in a gap region, the expression of a xy becomes simple. The Coleman-Hill theorems tells 
us that in this region only the lowest order diagram, i.e. the diagram of Fig. 4 with the bare propagator and the bare 
vertex, does contribute to a xy . Thus, the full propagator S can be replaced by the bare propagator Sq. The integral 

24rr 2 



J tr(dSo Sq) 3 gives a integer value and in fact counts the number of Landau bands below the Fermi energy. 



Thus, the Hall conductance is proved to be a integer times e jh in a gap region 



IV. PERIODIC POTENTIAL AND DEFECT 



In this section, we apply the formalism developed in the previous section to a system subjected to a periodic 
potential. Generally, a periodic potential lifts the degeneracy of the Landau level and the spectrum consists of 
Landau bands. As we mentioned in the introduction, Bloch electron in a strong magnetic field is of great interest in 
some contexts. In what follows, we assume that the potential lattice is formed by two linear-independent basis vectors 
with integer coefficients. Such a lattice is called a regular lattice. A periodic potential without this property, e.g. 
a honeycomb lattice potential, is not considered in this paper. A periodic potential with this property is generally 
written as 

V(x) = £ «(x + aN x wl ot + aN v wl ot ), (4.1) 

where wf f ,w^ of are basis vectors of the potential lattice. As is well-known, a periodic potential problem in a magnetic 
field is very sensitive to the flux penetrates a unit cell of the potential lattice. The unit cell means a parallelogram 
spanned by basis vectors of the potential lattice. If the flux $ is a rational multiple of the flux quantum $ , that is, 
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i=^=Im(<<)M ot = ^ (4-2) 

each Landau band splits into q subbands. The magnetic Brillouin zone (MBZ) is one q-th of the Brillouin zone. In 
addition, it was proven that the spectrum in the MBZ is p-fold degenerateE3. We re-prove this degeneracy in a simple 
manner using the previously devel oped formalism. 

When the flux is given by Eq. (4.2), it is convenient to select basis vectors of the von Neumann lattice as u) x = 
uj pot /q,uj y = pu> pot . The moduli of the von Neumann lattice becomes t = T pot /pq, where the moduli of the potential 
lattice is defined by r pot = —uj pot / ui pot . In this case, the potential energy term in the second quantized form becomes 



v 

p-1 q-l 



H *>* = ^EEE / d Mx) / ^b\{v)b v (p x -2* 



x{l,p x - 2tt(^ + -),p y + 2tt-|x = 0>(x = 0\l',p x - 2tt(- + - + -),p y + 2tt-). (4.3) 
a p a a p q a 

The eigenvalue equation reads 

(E - £,tylr(p) = - E E E / d 2 xv( X )(l, Px - 2n(V + - + -), Py + 2tt-|0) 

x(0\l',p x - 2tt(^ + - + -),p y + 27r-)Viv(p), (4.4) 
a p q a 

where p is a momentum in MBZ, that is, \p x \ < 7r/g, \p y \ < ir. The function V'/r(p) is given by the following form: 

7* 

^r(p) = in{Px - 2n-,p y ). (4.5) 

q 

The equation is solved by diagonalizing a Lq x Lq matrix, where L is the number of Landau levels. Therefore, each 
Landau band splits into q subbands generally. It is easy to see that the spectrum of the eigenvalue equation is invariant 
under the translations p x — > p x + 2irm/p and p x — > p x + 27rn/g, where m, n are integers. Since p, q are relatively prime 
integers, the above symmetry reads 

Ti 

E(p) = E{p x + 2Tr—,p y ). (4.6) 
pq 

Thus, it is proven that the spectrum in MBZ is p-fold degenerate. 

Next, we consider additional defects to a periodic potential. To extract properties in such a system, let us suppose 
a defect in a periodic short-range potential. That is, the potential is given by 

V(x) = E ^>a 2 <5(x - aN x w p x ot - aN y w pot ) + ga 2 6(x - aM x w p x ot - aM y ^ p y ot ), (4.7) 
Nez 

for integers M x ,M y . The second term breaks the periodicity of the potential. If we neglect the second term, the 
eigenvalue equation becomes 

E— P, t (p)A'(p))rr'#r'(p) = (E ~ ^)^r(p), (4.8) 

T?> q 

where p x q matrix Di is given by 

(D l (p)) sr =a(0\l )Px -2Tr(- + -) >Py ) O» = 0,..,p-l;r=0,..,g-l). (4.9) 

p q 

As was discussed in our previous papei0, the spectrum in the LLL consists of flat bands and Hofstadter-type bands. 
If we take account of higher Landau levels and the Landau level mixing, flat bands still exist at original Landau levels 
and also Hofstadter-type bands exist. The wave function of the flat band at E = Ei is given by tpk = ^«Ker(Dj). 







FIG. 6. The spectrum of the periodic short-range potential problem. The lowest two Landau bands are drawn in this figure. 
The lowest five Landau bands are taken into account. 



In Fig. 6 the spectrum for the square lattice potential is shown. Here, we incorporate lowest five Landau levels and 
Vo/frLu c = 0.3 as an example. We observe that a self-similar pattern and large gaps above flat bands exist in this figure. 
Hofstadter-type bands above flat bands tend to a set of bound states as t becomes infinity. Detailed explanation of 
the spectrum and the Hall conductance will be given in our forthcoming article. 



Incorporating the defect leads to the following additional term in the L.H.S. of Eq. (4 



I'm' 



^ 'mbz (27r) 2 



^p ; t (p)) rs (A'(q)) s ™'e l9Ma = (fe ^ )+lS(9 «- p « ) #m'(q) I (4-10) 



where Mod(M y ,p) = s and (M y — s)/p = S. Apparently, flat bands are still flat even if the defect exists. Furthermore, 
in each subband gap of the periodic short-range potential problem a bound state appears rearranging eigenstates of 
Hofstadter-type bands. The equation of bound state energies is given by 



V f #p |E ir (A(p)) s ,^/( P )^ M -^ + ^| 2 . 
9Z ?Jmbz(^) 2 E-E A (p) ' 1 ' > 



where ipi r A (p) and Ea(p) are the eigenfunction and the eigenvalue of Eq. (4.8). The R.H.S of the above equation 
generally becomes infinite when E approaches to the upper edge of the Landau subband. Also, it becomes minus 
infinite when E approaches to the lower edge. Thus, there exists a solution in each subband gap. The solution 
corresponds to a bound state trapped at the defect. If many defects exist in the periodic potential, many bound 
states appear in each the subband gap. As the number of defects increases, subband gaps tend to be filled with bound 
states and the number of extended states decreases. At last states which correspond to Landau subbands and bound 
states are shuffled each other and turn into a set of bound states (localized states) which has a finite density of state. 

V. QUANTIZED HALL CONDUCTANCE OF BLOCH ELECTRONS 

In this section, we generalize results given in section III to a system subjected to a periodic potential. In what 
follows, we assume the periodic potential has t — q/p. Thus, the each Landau band splits into q subbands generally. In 
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this case we can not directly use the Ward-Takahashi identity and thus the topological formula of the Hall conductance 
discussed in section III. This is because only the reduced momentum, i.e. the momentum in MBZ is conserved in this 
case. Therefore, we must generalize the Ward-Takahashi identity to the system subjected to a periodic potential first. 

Let us suppose that the eigenvalue and the eigenfunction of the one-body Hamiltonian under the periodic potential 
is given by Ea(j>), ipi r A (p). Here, p is a momentum in the MBZ, A is the suffix to specify the eigenstate, I is the 
Landau level index and r — 0,1, ...,q — 1. The electron propagator in the current basis is written as 

~ T t' 

Su>(p ,Px - 2ir-,p y ;p' Q ,p' x - 2-K—,p' y ) 

= (2^) 3 S(p-p')U lh (p x - 2n-,p y yj hr A (p) L^, 2r ,^*( p )tft - 2n-, Py ) 

q Pa - &a\P) 1 

= {2^f8{p~p')S [lrWr , ) {p). (5.1) 

Thus the momentum in the MBZ is conserved in the propagator. The propagator Sn r \n/ r /\(p) can be regarded as a 
matrix whose indices run over allowed set of (I, r). Owing to the orthogonality and the completeness of ip, its inverse 
is given by 

S-i = u^P^lA^ut. (5.2) 

i 

The one-particle irreducible vertex in the current basis is the same as the one in the free theory, because effects of the 
periodic potential are absorbed in the propagator. Thus, for an infinitesimally small momentum transfer, the vertex 
becomes 

~ T T 1 

^U>(Po,Px ~ 2ir-,p y ;p' ,p' x - 2ir—,p' y ;q) 
= -i8 rr , (2ir) 3 S(p + Q-p') x 

= (27r) 3 S(p + Q-p')f1 lrKl , rl) (p,p + Q). (5.3) 

Between S and T, the Ward-Takahashi identity is satisfied in a generalized form. The Ward-Takahashi identity in 
this case is simply given by 

The current correlation function for infinitesimally small momentum transfers q, q' is given by 

^{q,q') = \{2nY5(q-q') f^f -£?-Tr\r»(p,p + Q)S(p + Q)f" (p + Q ,p)S(p)} 
a J 2tt J mbz (2ny 

= {2^f5{ q -q')^{q). (5.5) 

Here, the trace is taken over indices (I, r). In the expression of the Hall conductance Eq. ( 0.15 ), the term including 
d p T p or d p Y v are vanished when the totally antisymmetric part is taken. The term including d p S survives. Therefore, 
using the Ward-Takahashi identity, the Hall conductance is written as a topologically invariant expression of the 
propagator S: 



1 



h 24tt 2 



TridS^S) 3 . (5.6) 

MBZxS 1 



Here, the p integral is taken over MBZ. Since the integral / Tr(dS~ 1 S) 3 gives the winding number of the 
propagator which is a integer value under general assumptions, the Hall conductance is proved to be a integer times 
e 2 /h whereas the filling factor is a fraction. This formula of the Hall conductance is a generalization of the topological 
formula Eq. ( |3.16| ) to the system subjected to a periodic potential. 
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Next, we consider the relation between the winding number and Chern number. Thouless et al. showed that the 
Hall conductance of Bloch electrons in a magnetic field is given by the Chern number of the Bloch function. This 
formula of the Hall conductance leads to a surprising result in a cosine potential. That is, in the gap region of 
Hofstadter bands, the Hall conductance becomes a integer times e 2 /h in spite that the system has a fractional filling 
factor in the gap. Furthermore, the Hall conductance is changed drastically as the Fermi energy is increased across 
subband gaps. This shows each subband carries a very ilyge mobility. These features of the Hall conductance in a 
cosine potential are also obtained by the Streda formulanJ. In contrast to the cosine potential the Hall conductance 
in a periodic short-range potential is not known. Taking account of recent antidot experiments, it is also important 
to study the Hall conductance in the periodic short-range potential. 

To solve these problems, we rewrite the generalized topological formula Eq,-( j5.6| ) in a more convenient form. In 
computation of Eq. ( |5.6| ), there is a useful relation due to Polyakov-WiegmannEa 

I (ST) = 1(g) + I(T) + 3 J TrdiS^dSdTT- 1 ), (5.7) 

where I(S) is defined by J Tr(SdS~ 1 ) 3 . To utilize the relation, it is convenient to regard ip, U, = i/(po — EjCj as 
matrices of order Lq. We regard that Uw (p x —2nr/q, p y ) is diagonal with respect to the index r and that is diagonal 
with respect to the ind ex A . Therefore, the propagator S is regarded as a product of three matrices, ip u = Uip, 



and (ip u y ■ Using Eq. (5/7), the generalized topological formula leads to 

'--ts £ (L*" V) - 2 "i)' (5 ' 8) 

A(E A <E F ) H/ 

This formula is similar to that of Thouless et aM and Kohmotoil except that there is the second term proportional to 
the filling factor. The second term comes from the matrix U(p). However, this term at last cancels with the boundacy 
contribution of the first term. For the first term, we give a topological argument which is similar to the one given in0. 
The boundary condition in BZ implies that the phase of the eigenfuction varies by 2ir as one follows the boundary of 
BZ counterclockwise. Since BZ is divided into q MBZs, this leads to a constraint given by 

9-1 . 

V <L darg(^/) = 2tt. (5.9) 

r=0 J dMBZ 

Under this constraint the configuration which is singular in ip A dip A at some p may appears for a band A. We call 
this the vortex. If the vortex appears, ip^dip is ill-defined in the entire MBZ. Let us assume there is a vortex at 
p = p c . We remark that there is a gauge degree of freedom in the eigenfunction. That is, the overall phase factor 
of the eigenfunction ip A (p) is n °t determined. As is well known in the description of monopole, we can define ^ dip 
in the entire MBZ using the gauge degree of freedom. In order to define it, we divide MBZ into two regions. One 
region (MBZi) does not include p c and has a boundary of MBZ. The other region (MBZ2) includes p c . The boundary 
between two regions is a closed line around p c . In MBZ2 using the gauge degree of freedom, we choose a different 
phase assignment which has no vortex. It is not necessary to satisfy the boundary condition when we extrapolate the 
wave function to MBZi. Thus, ip^dip is completely defined in the entire MBZ. However, there is a phase mismatch at 
the boundary between two regions. That is, the eigenfunction tp^ 1 ' in MBZi and the eigenfunction in MBZ2 are 
related as = e t6 ip^ at the boundary. For the integral / MBZ . d(ip^dip)(i — 1,2), we can use the Stokes theorem 
and obtain 

d(ipUip)= [ ipUip + i [ dO, (5.10) 

MBZ J 9MBZ Jc 

where C is the closed line around p c . The first term in R.H.S. of the above equation is equal to 2iri/q due to the 
normalization and the boundary condition of the eigenfunction ip. The second term gives 2iri times a integer. This 
integer is called the Chern number. Therefore, starting from the generalized topological formula we obtain the Chern 
number for the Hall conductance: 

A(E A <E F ) JO 

To compare our result with the one obtained by Thouless et aZ.§, we consider the correspondence between the Bloch 
function itk and the eigenfunction tpi r in the von Neumann lattice basis. In terms of the Bloch function the energy 
eigenstate can be written as 
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«k(x)e lkx . 



(5.12) 



On the other hand, in the von Neumann lattice basis the energy eigenstate can be written as 

1 x r 

—^22^lr(p)(^\l,Px-^-,Py)- (5.13) 

yQ " 



Lr 



Both functions satisfy the same boundary condition (except for a constant phase) and normalization in the magnetic 
unit cell (MUC) if the momentu m k i s replaced by the momentu m p w ith the rule ak • wg°* = qp Xl pak ■ w^ * = p y 



in a t = q/p problem. Thus Eq. (5.12) can be identified with Eq. (5.13). The Chern number obtained by Thouless et 
al. is given by 

n c = ^-.\ dA, (5.14) 

where A is the connection 1-form written as 

A(k) = / d 2 xuk(x)d k u k (x). (5.15) 
Jmuc 



Substituting Eq. (5.13) into Eq. (5.14), we obtain 

1 f I 



nc 



dtyUip) - -. (5.16) 



27r« Jmbz <7 



This Chern number is nothing but the one in Eq. (5.8). Therefore, the winding number of the propagator is reduced 
to the Chern number of Thouless et. al. if interactions are neglected. 

So far, we neglected interactions and consider only a periodic potential. However, in an interacting system sub- 
jected to the periodic potential the winding number formula of the Hall conductance can be applied in the following 
way. In this system we can treat the perio dic potential nonperturbatively as discussed above and treat interactions 



perturbatively. Thus, the propagator Eq. (|5.l| ) replaces the "bare" propagator. Since the Ward-Takahashi identity 



and related theorems still hold, results obtained in Sec. 3 can be also applied in th is s ystem. Therefore, the Hall 



conductance in a gap region is written by the winding number of the propagator Eq. (5.1). Furthermore, the winding 
number is the same as the Chern number obtained by Thouless et al. as discussed above. 

Finally we calculate the Hall conductance in the system of a periodic short-range potential. To calculate the Chern 
number for a general periodic potential is a subtle problem. However, we can calculate the Chern number for some 
analytically solved potentials. Let us first consider the potential with t — 1/p. This is the case that each Landau 
band does not split into subbands. In this case, the p-dependence of ipi(p) can be set as /3(p)/|/?(p)| which satisfies 
the boundary condition. Apparently, this has a vortex at p = 0. Thus, we must choose another gauge for the wave 
function around p = 0. A simple choice is 1, that is, there is no p-dependence around p = 0. Since § c darg/3(p) = 2tt, 
the Chern number of each Landau band is equal to 1. Thus, the Hall conductance is e 2 /h times the number of Landau 
bands below the Fermi energy as expected. 

Next we consider a t = 2 case. A periodic potential with t = 2 can be solved analytically if the number of 
truncated Landau levels is one or two. Let us consider the periodic short-range potential in the lowest Landau level 
approximation. In this case the spectrum consists of a flat band and a lifted band. The wave function of the flat band 
satisfies the equation 

P*(p)Mp) + P*(p* - k,p v )Mp) = o. (5.17) 

Here, ipk(p) = ip{p x — kn,p y ) and MBZ is the region \p x \ < tt/2, \p y \ < n. The following choice of the phase is 
consistent with the boundary condition of y-direction: 



^o(p) \ = 1 ( \/3( Px - ir,p y )\e la ^P»+ l P* 

MP) ) ^a(p)+a( Px -n,p y ) \ -|/3(p)|e*«^-- 



(5.18) 



At p = the upper component of the above eigenfunction has a vortex. Therefore, we must choose another patch 
around p = 0. A possible choice of the eigenfunction around p = is given by 



MP) \ _ 1 ( P*{Px ~ TT,p y ) 

MP) J y/a(p) + a(j> x -n,p y ) V -0*(P) 



(5.19) 
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Thus there is a phase mismatch = — arg/3(p) — axg/3(p x — n,p y ) at the boundary. Its integral along the boundary 
is — 2ir. Thus, we conclude that the Hall conductance of the flat band is e 2 /h. The eigenvalue equation for the lifted 
band becomes 

p(p x - ir, Py )Mp) ~ P(p)Mp) = 0. (5.20) 
Its solution which is consistent with the boundary condition is 



Mp) \ = 1 f /?(P) 

V'i(p) J v/a(p) + a(p x - ir,p y ) V P(P X ~ K'Py) 



(5.21) 



This solution has no vortex. Thus, the Hall conductance of the lifted band is 0. We can see that the Wannier state 
constructed from the flat band is relatively spread than the one from the lifted band. 

In the case that t is an integer larger than 2, the spectrum also consists of a flat band and a lifted band in the 
lowest Landau level. The eigenfunction of the lifted band 

MP) = , 2t=MM (5.22) 



satisfies the boundary condition and has no vortex in the MBZ. Thus, the Hall conductance of the lifted band is 0. 
This result is reasonable because the lifted band becomes a set of bound states and thus has no mobility as t becomes 
infinity. Since it is proved that the Hall conductance of the Landau band is e 2 /h, the Hall conductance of the flat 
band is e 2 /h. These results are irrespective to the lattice structure of the periodic short-range potential. 



VI. SUMMARY AND DISCUSSION 



We constructed three sets of one-particle states based on the von Neumann lattice in a gauge- independent manner. 
The first one is the set of coherent sates which are localized on the sites of the von Neumann lattice and not orthogonal. 
The second one is the set of momentum states which are extended, orthogonal and normalized with 5-function. The 
third one is the set of Wannier states which is localized on sites of the von Neumann lattice and is orthonormal. Using 
the second one we developed a field-theoretical formalism of a quantum Hall system and showed the Ward-Takahashi 
identity and the topological formula of the Hall conductance. This formalism was applied to a system subjected to a 
periodic potential. The p-fold degeneracy of the spectrum in a t = q/p problem is easily proved using this formalism. 
The generalized topological formula for the Hall conductance of Bloch electron was obtained. The formula is written 
as a winding number of the propagator. Relation between the winding number and the Chern number was clarified. 
As an example we calculated the quantized Hall conductance in a periodic short-range potential. 

In this paper we study the system subjected to a periodic potential using the set of momentum states. We believe 
that cuir formalism is also useful to study Coulomb interaction in the quantum liquid phase 2 !] and the Wigner crystal 
phasazj of the quantum Hall system. In the Wigner crystal phase, a triangular lattice potential is generated from 
the mean field of the charge density. In this case, the filling factor v is related to t with v — 1/t. If the mean field 
contributes to the short-range potential with positive coefficient, the Fermi energy lies in the large energy gap at 
v = 1 — 1/t. That is, the flat band of the lowest Landau level in Fig. 6 is completely filled. Then, v = 1/2 state can 
be a self-consistent solution which has a large energy gap in contrast to the quantum liquid phase in which gap closes 
at v = 1/2. If the mean field contributes to the short-range potential with negative coefficient, the Fermi energy lies 
in the large energy gap at v = 1/t. That is, the lifted band is completely filled. Then, all states of v < 1 can be a 
self-consistent solution which has a large energy gap. Since the Hall conductance of the lifted band is zero, the ground 
state in this approximation is close to the crystal. 

Although the short-range potential is not a good approximation for the Coulomb interaction, the mean field treat- 
ment starting from the short-range potential may give some insights to the Wigner crystal phase of the quantum Hall 
system. 
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APPENDIX A: 



The properties of functions used in this paper are given here. A theta function ■&i(z\t) is defined by 

oo 

&i{z\t) = 2 ^(-l)" e 7r "("+5) 2 sin[7rz(2n + 1)]. 



(Al) 



Under the modular transformationEZI, /3(p|r) defined in Eq. ( 2.15 ) transforms as follows: 



ct + a y \ct + a] 

x (5{dp x + bp y , qp x + apy\r), 



(A2) 



° b d )tSL(2,Z), 



where r\ is an eighth-root of unity. a(p) defined in Eq. ( 2.14 ) obeys the periodic boundary condition. G = In a is a 
Green function of the Laplacian operator on a torusES, that is, 



AG(z) = 2n5 {2) {z) 



47T 

Imr ' 



(A3) 



where 



6^(z) = ^2s^(z + N x + rN v ), 



N 



z = (Px + T~p y )/2ir, and A = 4d z d z . As mentioned in Sec. II, the isolated singularity of the delta-function at p = 
in Eq. (A.3) is not involved in our theory. The factor (fi\e~ tk '^\fi') which appears in p(k) is written by 



(/^- ^fe • ^ l/^) = ^/55(- £^ w 3£l )' 



xi{ n) (¥), 



(A4) 



where Lj is the associated Laguerre polynomial. 
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